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Nonlinear systems can be probed by driving them with two or more pure tones while measuring the intermod- 
ulation products of the drive tones in the response. We describe a digital lock-in analyzer which is designed 
explicitly for this purpose. The analyzer is implemented on a field-programmable gate array, providing speed 
in analysis, real-time feedback and stability in operation. The use of the analyzer is demonstrated for Inter- 
modulation Atomic Force Microscopy. A generalization of the intermodulation spectral technique to arbitrary 
drive waveforms is discussed. 



Nonlinearity is essential for the function of digital cir- 
cuits and amplifiers, and yet it is often considered the 
bane of electronic systems, causing undesirable frequency 
mixing known as intermodulation distortion (IMD). In 
some cases however, IMD can be useful for signal 
conditioning, for example in the phase-sensitive non- 
degenerate parametric amplifier, where a strong pump 
tone is used to deliver power gain to a weaker signal 
toneii 2 -. IMD can also be considered as useful if one 
is trying to gain understanding about the nonlinear- 
ity itself, for example in Intermodulation Atomic Force 
Microscopy 3 . Whether one wants to use or avoid non- 
linearity, a complete understanding of the nonlinearity 
is often desired, and to this end we developed a digital 
lockin analyzer for probing the real-time intermodulation 
response of a nonlinear two-port system. 

The modeling and parametrization of nonlinearity 
is the object of current microwave instrumentation 
development, where there exists an engineering need for 
accurate models of the large amplitude behavior of mi- 
crowave devices. The vector network analyzer, the mi- 
crowave version of the lockin, has recently been adapted 
for this purpose to measure the so-called X-parameters. 
This measurement is based on the poly-harmonic dis- 
tortion models, which assumes the existence of a strong 
drive tone, and approximates the mixing of weaker sig- 
nals which are harmonics of the drive tone with a linear 
mapping in frequency space. The lockin analyzer dis- 
cussed here is designed to measure intermodulation, or 
the frequency mixing of two or more drive tones applied 
to a nonlinear system. 



LOCKIN MEASUREMENT, HARMONICS AND 
INTERMODULATION 

The lockin measurement technique is defined by its 
use of a reference oscillation for making a narrow-band 
measurement of some physical process^. The digital in- 
termodulation lockin analyzer (ImLA) described here is 
capable of performing multiple lockin measurements si- 
multaneously and for such measurements, clock synchro- 
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FIG. 1. The frequency spectrum of two pure tones is mea- 
sured by sampling and discrete Fourier transform. The large 
background near the peaks (red curve) arises from the re- 
moval of synchronization between the clock in the waveform 
synthesizer and the sampling clock. 



nization is very important. Figure [T] compares two mea- 
surements, where two sine waveforms were added in a 
summing amplifier, sampled and Fourier transformed. 
We see that the removal of the synchronization between 
the synthesis clock and the sampling clock causes the ap- 
pearance of Fourier leakage. This leakage dramatically 
raises the measurement floor in the vicinity of the drive 
frequencies, below which real signals with good signal-to- 
noise ratio could be hiding. Synchronization establishes 
one master clock from which a common reference signal 
is derived, and is thus a key aspect of digital lockin mea- 
surements. 

The lockin method is often used to probe the response 
of a two-port system to a pure drive tone. The drive ap- 
plied at one port, Vd cos(wd£), is used as a reference for 
the measurement of two channels of response, V(t), at 
the other port. The lockin measurement is then realized 
by multiplying two unity-amplitude copies of the drive 
signal, one being phase-shifted by 7r/2, and integrating 
each for a time T. The result will be the Fourier co- 
sine and sine components of the response in a frequency 
bandwidth 1/T, centered at the reference frequency wd, 
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In the digital world we can realize this measurement 
by sampling the response at discrete times and calcu- 
lating a Fourier sum. The cosine and sine quadratures 
V x and V y can be alternatively represented as amplitude 

\V\ = yVj? + Vy and phase 9, where tand = V y /V x . 

The reference signal is taken to have zero phase, and the 
response phase 9 is the phase with respect to this zero. 

This approach to measuring signals is natural for linear 
systems driven by a pure tone, where the response will 
only be at the drive frequency. For nonlinear systems the 
response to one pure drive tone will include harmonics, 
or response at integer multiples of the drive frequency. 
The lockin amplifier can be easily extended to include 
harmonic analysis by simply forming two new measure- 
ment channels for each harmonic frequency desired. In 
this manner, the phase of the harmonics with respect to 
the reference signal is unambiguously defined. 

When a nonlinear system is driven with two pure tones 
at ui\ and W2, the response contains not only harmonics of 
the drive tones, but also intermodulation products (mix- 
ing products) of the drive tones, occurring at frequencies 
which are integer linear combinations of the two drive 
frequencies, loimp — fci^i + ^2^2 where ki,k 2 G Z. For 
the intermodulation response, which does not occur at 
a drive frequency or a harmonic, there is no reference 
signal and an unambiguous definition of phase cannot 
be found. We can however measure the phase of inter- 
modulation products by restricting our measurement so 
that all drive tones occur at frequencies which are integer 
multiples of some base frequency, which we call Auj. In 
this case all intermodulation products occur at integer 
multiples of Aw, and Auj can be used as the reference 
frequency in the lockin calculation. 



INTERMODULATION LOCKIN ANALYZER 

We used a field-programmable gate array with one 
channel of A/D input (12 bit) and two channels of D/A 
output (16 bit) to implement an intermodulation lockin 
measurement, where all relevant frequencies are integer 
multiples of Auj. Schematically depicted in hg.[5J the In- 
termodulation Lockin Analyzer (ImLA) synthesizes two 
user-dehncd drive tones and calculates both quadratures 
of the intermodulation response at 32 user-defined in- 
termodulation frequencies. The user communicates with 
the ImLA via an ethernet connection so the instrument 
can easily be integrated into measurement systems using 
a wide variety of software platforms (e.g. LabView®, 
MATLAB®, Python®, etc.). 
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FIG. 2. A schematic showing the functionality of the digi- 
tal Intermodulation Lockin Analyzer (ImLA). The frequencies 
(fi), drive amplitudes (Ai, A2), feedback gain (P) and bias 
(Vt>) can all be controlled via the ethernet interface. For each 
of the 32 frequencies both the inphase (J) and the quadrature 
(Q) component is generated and used in a lockin calculation. 



In its present form the ImLA samples at 61.4 MSa/s. 
64 Fourier sums are calculated in parallel by accumulat- 
ing every sample during a period of time T, which is the 
inverse of the user-defined measurement bandwidth 1/T. 
The quadrature response at 32 user-defined frequencies 
is thus acquired and sent over ethernet to a computer. 
For one of the 32 frequencies the amplitude and phase 
are calculated by the CORDIC algorithm! in the FPGA, 
so that these signals can be used for real-time feedback. 
The update rate of the feedback can be selected up to 
1024/T. This feature is useful in AFM systems where 
fast feedback is needed to control the scanning probe's 
height above the surface. 

The Intermodulation lockin can also be run in a time- 
domain mode, where the response is continuously down- 
sampled to 3.9 MSa/s and streamed over ethernet with- 
out gaps in the data. This mode is useful for detailed 
study of the entire response spectrum. For synchroniza- 
tion with other measurement instruments, the ImLA has 
a 10 MHz clock OUT signal and an external trigger in- 
put which will reset the calculation of the Fourier sums, 
effectively defining the start of measurement. Many of 
the ImLAs specific features and performance specifica- 
tions can be modified by reconfiguration of the firmware 
on the FPGA. 

The ImLA's ability to measure nonlinearity in a De- 
vice Under Test (DUT) is limited by its own internal 
nonlinearity We tested this limit by sending a maxi- 
mum amplitude drive signal directly from the OUT port 
to the IN port, to measure an internal Total Harmonic 
Distortion (THD) and IMD, which were both less than 
—75 dB. Under these test conditions we also find that 
the dynamic reserve at 1 kHz measurement bandwidth is 
80 dB. 
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INTERMODULATION ATOMIC FORCE MICROSCOPY 

Figure [3] demonstrates the use of the ImLA in Inter- 
modulation Atomic Force Microscopy 3 where IMD in a 
mechanical oscillator is measured. Typically the first 
flexural eigenmode of the cantilever is excited, which of- 
ten has a sharp resonance around 350 kHz with quality 
factor Q ~ 600 in air. The ImLA is used to drive the can- 
tilever with two tones separated by ~ 500 Hz, centered 
around resonance. The cantilever response is detected by 
optical means in the AFM, and the output of a photo- 
diode circuit is analyzed by the ImLA. 

The free cantilever oscillating well above the surface 
is a linear system and the response to two drive tones 
closely spaced in frequency will be a beating waveform in 
the time domain (fig. [3] shows one beat). When the oscil- 
lating cantilever approaches an SiC>2 surface, the sharp 
tip on the end of the cantilever experiences a nonlinear 
tip-surface force, which creates intermodulation of the 
drive tones. Figure [3] shows two frames from a movie 
of an approach-retract measurement, which can be seen 
online-. In these measurements the ImLA was used in 
time-domain mode and the spectra were calculated by 
Fourier transform after the measurement. 

When the AFM is scanning over a surface, speed in 
analysis is required, and the ImLA is used in lockin mode, 
where the intermodulation response is captured only at 
the peaks in fig. [3] (lower right panel) . While scanning, 
multiple images can be viewed by plotting the amplitude 
or phase at the various intermodulation frequencies as 
a function of tip position™ From the intermodulation 
spectrum measured at each pixel of the image, one can 
reconstruct the tip-surface interaction forced, and reveal 
the nano-scale mechanical properties of the surface. This 
nonlinear interaction force is the information which can 
be measured with dynamic AFM, and one can view the 
IMD measurement not only as a sensitive probe of this 
nonlinear interaction, but also as an efficient means of 
compressing information in the sampled data. By lock- 
ing in on the intermodulation frequencies we extract for 
storage only the essential amplitudes and phases in the 
response, disregarding only noise. 

MULTI-FREQUENCY INTERMODULATION ANALYSIS 

Thus far we have considered a drive consisting of two 
pure tones and a measurement of the intermodulation 
products of these two tones. A natural extension of 
this technique would be to include arbitrary drive wave- 
forms. For example, we could synthesize a drive wavelet 
other than a simple beating waveform, such as a chirplet 
formed by combining several frequencies in a band 11 . 
The drive waveform should be built from a superposi- 
tion of pure tones where the frequencies are taken from a 
set of frequencies formed by the integer multiples of Aw. 
A lockin method will thus be possible since intermodula- 
tion response will occur only at integer multiples of Aw. 
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FIG. 3. Time domain (upper panels) and frequency domain 
(lower panels) of free cantilever (left panels) and an engaged 
cantilever (right panels) in Intermodulation AFM (enhanced 
online 8 ). The time domain shows the sampled response over 
only one beat of the drive wavelet. The frequency domain 
shows the FFT of the sampled data over 10 consecutive drive 
beats. The small measurement bandwidth used in the FFT 
allows one to observe how the noise level in between the in- 
termodulation peaks changes while approaching the surface. 

In this case, the analysis of the response is identical to 
that described here, and the architecture of the ImLA 
is well suited for such an extension. However, for mul- 
tiple drive tones, the measured intermodulation would 
include all possible mixing products of the many drive 
frequencies. Nevertheless, the theoretical approach based 
on polynomial representation of the nonlinearity to cre- 
ate a linear mapping of the intermodulation products in 
frequency spaced should be able to handle this compli- 
cation. It is possible that suitably chosen drive wavelets 
could be advantageous for the intermodulation spectral 
technique, depending on the application and the infor- 
mation one wishes to extract from the intermodulation 
measurement. 

In summary, we described a digital incarnation of a 
lockin measurement instrument specifically designed to 
measure the intermodulation of pure tones by nonlinear- 
ity. This Intermodulation Lockin Analyzer (ImLA) can 
find use in many measurement applications where one 
wishes to make a detailed model of the nonlinearity of 
a device under test. The intermodulation spectral tech- 
nique is particularly useful in applications where the per- 
turbation of a high quality factor resonance is employed 
for sensing. We demonstrated one such an application in 
the field of dynamic Atomic Force Microscopy, however 
numerous sensors based on this measurement principle 
could use the intermodulation lockin to gather more in- 
formation than is otherwise possible with standard lockin 
techniques based on a single drive frequency. 
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